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ABSTRACT 

We present new measurements of the Rossiter-McLaughlin (RM) effect for three WASP 
planetary systems, WASP-16, WASP-25 and WASP-31, from a combined analysis of 
their complete sets of photometric and spectroscopic data. We find a low amplitude 
RM effect for WASP-16 (Teff = 5700 ± 150 K), suggesting that the star is a slow 
rotator and thus of an advanced age, and obtain a projected alignment angle of A = 
-4.2°+ 1 1 3 1 9 °. For WASP-25 (Teff = 5750±100K) we detect a projected spin-orbit angle 
of A = 14.6° ± 6.7. WASP-31 (Teff = 6300 ± 100 K) is found to be well-aligned, with a 
projected spin-orbit angle of A = 2.8° ±3.1. A circular orbit is consistent with the data 
for all three systems, in agreement with their respective discovery papers. We consider 
the results for these systems in the context of the ensemble of RM measurements made 
to date. We find that whilst WASP-16 fits the hypothesis of iWinn et al.M2010|) that 



'cool' stars (Teff < 6250 K) are preferentially aligned, WASP-31 has little impact on 
the proposed trend. We bring the total distribution of the true spin-orbit alignment 
angle, "0, up to date, noting that r ecent results have improved the agreement with the 
theory of Fabrycky & Tremaine (|2007[) at mid-range angles. We also suggest a new 



test for judging misalignment using the Bayesian Information Criterion, according to 
which WASP-25 b's orbit should be considered to be aligned. 

Key words: planetary systems - stars:individual: WASP-16 - stars:individual:WASP- 
25 - stars:individual:WASP-31 - techniques:radial velocities 



1 INTRODUCTION 

As the number of transiting "hot Jupiters" known to astron- 
omy has grown, there has been a gradually increasing push 
towards fully categorising their physical and orbital proper- 



* based on observations made using the CORALIE high reso- 
lution echelle spectrograph mounted on the 1.2 m Euler Swiss 
Telescope and the HARPS high resolution echelle spectrograph 
mounted on the ESO 3.6 m (under proposals 084.C-0185 & 085. C- 
0393), both at the ESO La Silla observatory. 
f E-mail:djab@st-andrews. ac.uk 



ties. It is widely presumed that close-in gas giants do not 
form at the locations in which we observe them, and there 
are competing theories to describe the process that leads 
them to their observable orbits. 

Migration induced by a protoplanetary disc provides 
one means by which such a situation can be explained 



(Lin, Bodenheimer & Richardson 19961. Since such discs 



are generally aligned with the host star owing to angular 
momentum conservation, we would expect that disc migra- 
tion would preferentially produce well-aligned hot Jupiter 
systems. Some misaligned planets would not be out of place 
under this mechanism, being the result of close planet-planet 
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encounters following migration, but we would expect the ma- 
jority of planets to exhibit spin-orbit alignment. 

The Kozai-Lidov mechanism | Kozai|[l962 Lidov|[l962 1 
is the basis of a competing theory for which evidence is 
mounting. The presence of a third, outer body in a plan- 
etary system can excite periodic oscillations in both the ec- 
centricity and inclination of a planetary orbit; inward migra- 
tion then follows, with tidal friction kicking in as the planet 
approaches its host, causing the orbit to shrink and circu- 



by Winn et al. (20101 (hereafter W10I, who speculated 



larise (Fabrycky & Tremaine 20071. The oscillating incli- 



nation that results from Kozai-Lidov interactions produces 
a continuum of inclinations once the orbits are stable, and 
thus we would expect the majority of hot Jupiters to exhibit 
misaligned orbits if the Kozai-Lidov mechanism operates. 

It is possible, to some extent, to distinguish between 
these competing theories through measurement of the spin- 
orbit alignment angles of hot Jupiter systems. Given the 
different angular distributions predicted by these theories, 
building up a significant number of spin-angle measurements 
is a useful means of determining which mechanism is act- 
ing. Unfortunately the true misalignment angle cannot be 
measured unless a spectroscopic measurement of v sin I is 
made, and the stellar rotation period is known. This yields 
an estimate of the inclination axis to the line-of-sight (e.g. 



Schlaufman ( 2010 1) . The situation is made more difficult by 
systematic uncertainties in vsinl measurements, and the 
sine function, which flattens as it approaches 90° and there- 
fore only yields useful measurements at low to intermedi- 
ate inclinations. We are thus currently limited to measuring 
he projected spin-orbit misalignment angle in the plane of 
the sky. This is generally measured through the Rossiter- 
McLaughlin (RM) effect ( |Rossiter||T924} |McLaughlin|[T924 1 
which is observable during transit. As the planet transits 
the approaching limb of the star its spectrum is red-shifted, 
and when it transits the receding limb its spectrum is blue- 
shifted. The precise form of the RM anomaly in the radial 
velocity (RV) curve gives the projected misalignment angle, 
A. 

The first observation of the RM effect for a transiting 
planet was made by |Queloz et al. ( 2000[) , and since then 



the number of measurements has increased significantly to a 
level such that it is possible to begin carrying out analysis of 



the ensemble of measurements. Fabrycky and Winn ( 2009 1 



investigated 11 systems with known values of A, deriving two 
theoretical distributions for i\>, the true misalignment angle, 
using different assumptions about the form of the distribu- 
tion. They suggested, based on an apparent dual population 
within their data set, that there might be two routes for 
planet migration, one producing mostly aligned planets and 
the other producing misaligned planets. 

One early indication of a pattern was that misaligned 
planets tended to be high mass and on eccentric orbits 
( jJohnson et al.|[2009| ). Subsequent observations have often 
countered this initial trend (for example HAT-P-7 ( Winn et 
al.|2009|[Narita et al.|2009| ), but high mass (M P > 4MJup) 
planets do appear to have a different obliquity distribution 
( Hebrard et al.|2011 I . Of the 6 planets in this category with 
measured misalignment angles four are misaligned, but none 
have |A| > 50°. More observations of high mass planets are 
needed before we can be certain that this is not merely an 
artefact of small-number statistics however. 

One of the more intriguing suggestions was put forward 



that the division into aligned and misaligned planets might 
be dependent on the effective temperature of the host star. 
Using a larger sample of 19 systems with known A, they 
found that the misaligned systems were preferentially hot- 
ter than the aligned examples, with a critical temperature of 
Teff ~ 6250 K dividing the two populations. One explana- 
tion put forward for this was the tidal realignment of planets 
around 'cool' stars, with the equivalent process around 'hot' 
stars being suppressed owing to their lack of a convective 
envelope. |W10| further conjecture that the current ip distri- 
bution could be completely explained by a migration mech- 
anism driven by a combination of Kozai-Lidov oscillations 
and planet-planet scattering, without the need to invoke disc 
migration. 



Triaud et al. (20101 (hereafter TlOl added 6 planets 



to the ensemble of known RM measurements. Calculating 
individual ip distributions for each planet based on the as- 
sumption that stellar rotation axes are randomly oriented 
on the sky, they produced a total distribution for the en- 
semble of planets, finding that it matched the theoretical 
distribution of Fabrycky & Tremaine (2007) for Kozai-Lidov 



mechanism dominated migration, further implying that disc 
migration might be superfluous to requirements for explain- 
ing the presence of hot Jupiters. 

Here we present measurements of the RM angle for 
three more planets from the Wide Angle Search for tran- 
siting Planets (WASP) ( |Pollacco etirpOOB) , WASP-16b, 
WASP-25 b and WASP-31 b, and investigate how they mod- 
ify the ensemble results and conclusions discussed above. In 
section[2]we give details of our observations, and in section[3] 
we discuss the analytical methods used to determine the 
misalignment angles. In sectionHlwe report on the results of 
our analysis for the individual systems. In section|5]we dis- 
cuss the implications of our results for previously observed 
trends. Finally, in section[6] we take another look at the 
question of alignment, presenting a new test for planetary 
orbit misalignment. 



2 OBSERVATIONS 

Radial velocity data for all three planetary systems were 
obtained using the CORALIE high precision echelle spec- 
trograph ( [Queloz et al.|2000b l mounted on the Swiss 1.2 m 
Euler telescope, and with the HARPS high precision echelle 
spectrograph (May or et al.|200"3| mounted on the 3.6 m ESO 
telescope at La Silla. Data from CORALIE were used pri- 
marily to constrain the presence of a long-term trend in ra- 
dial velocity that might be indicative of a third body in the 
system, whilst HARPS was used to monitor the radial ve- 
locity before, during and after a specific transit event. Two 
data points were obtained the night before the transit, and 
for at least one night following the transit; on the night of 
the transit observations were started 90 minutes prior to the 
predicted start of transit and continued until 90 minutes af- 
ter its predicted conclusion. 



2.1 WASP- 16 

WASP- 16 was observed using CORALIE between 2008 
March 10 and 2009 June 3, on an ad-hoc basis. One datum 
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Table 1. System parameters for the three WASP planetary systems for which we evaluate the Rossiter-McLaughlin effect. Parameters for 
WASP-16 were taken from|Lister et al.H2009) . Parameters for WASP-25 were taken from |Enoch et al.| < ]2011) . Parameters for WASP-31 
were taken from | Anderson et al.|||2011|. usi n I and macroturbulence values have been updated through spectroscopic analysis of the new 
HARPS data using the|Bruntt et al,|||2010|l calibration. 



Parameter 



Unit 



WASP-16 



WASP-25 



WASP-31 



M* 


M @ 


1 022+ ' 074 


ri# 


Re 


Q46+ 0057 

U.94D_ 052 


Teff 


K 


5700 ± 150 


v sin I 


kms" 1 


2.3 ±0.4 


macroturbulence 


krns" 1 


2.3 


M v 


MJup 


u.ooo_ 076 


Rp 


-RJup 


1 nns+ 0083 
i.uo»_ 0060 


P 


days 


3.11860 ±0.0( 


a 


AU 


0421+ 00010 
u.u^zi_ 0019 


e 




O(adopted) 


i 





OCT QQ+0.27 



1.00 ±0.03 
0.92 ±0.04 
5750 ± 100 
2.6 ±0.4 
2.4 
0.58 ±0.04 

, 99 +0.06 

-■-• zz -0.05 

3.764825 ± 0.000005 
0.0473 ± 0.0004 
O(adopted) 
88.0 ±0.5 



1.161 ±0.026 

1.241 ± 0.039 

6300 ± 100 

8.1 ±0.5 

1.2 

0.478 ± 0.030 

1.537 ±0.060 

3.405909 ± 0.000005 

0.04657 ± 0.00034 

O(adopted) 

84.54 ±0.27 



was also acquired on 2010 July 14 to retest the hypothesis 
of a long-term radial velocity trend. The transit observed 
with HARPS occurred on the night of 2010 March 21; 32 
data points were acquired over the duration of the night. 
This transit observation was affected by cloud cover, so an 
additional transit was observed on the night of 2011 May 
12, producing a further 28 RV measurements. Further mea- 
surements were made on the days surrounding this transit 



3 DATA ANALYSIS 



as well (see journal of observations, Tables Bl B2 and B3l 



Details of the photometric observations of WASP-16 are 
given in |Lister et aT] ( |2009[ ) . 



2.2 WASP-25 

HARPS observed the transit taking place on the night of 
2008 April 11. 44 observations were made that night, with 
additional data acquired on adjacent nights (see the jour- 
nal of observations, Tables B4 and B5 I. The system was ob- 



served using CORALIE between 2008 December 29 and 2009 
June 28, with observations made at irregular intervals be- 
tween these dates. 

|Enoch et al.| ( |2011[ ) describe the photometric observa- 
tions that were made of WASP-25. 



2.3 WASP-31 

WASP-31 was observed using CORALIE between 2009 Jan- 
uary 4 and 2010 May 18 during several short runs. HARPS 
was used to observe a full transit on the night of 2010 April 
15, with 17 data points obtained. Additional observations 
were made on adjacent nights (see the journal of observa- 
tions, Tables|B6]and[B7l. 

The photometric observations for WASP-31 are dis- 



Our analysis mirrors that of |T10[ using the complete set of 
photometric and spectroscopic data for the objects that we 
investigate in order to fully account for parameter correla- 
tions. We use an adapted version of the code described in 
|Collier Cameron et al.| ( |2007[ ), fitting models of the photo- 
metric transit, the Keplerian RV and the RM effect to the 
system data. The fit of our model is refined using a Markov 
Chain Monte Carlo (MCMC) technique to minimize the \ 2 
statistic, and to explore the parameter space using the jump 
parameters To (epoch of mid-transit), P (orbital period), 
W (transit width), 6 (impact parameter), 7 velocity, 7, K 
(RV semi-amplitude), Teff (stellar effective temperature), 
[Fe/H] (metallicity), -^/ecosoj, -^/esino;, \/v sin I cos A and 
\/v sin /sin A. We use a burn-in phase of 2000 steps, with 
burn-in judged to be complete when \ 2 becomes greater 
than the median of all previous values (Knutson et al. 2008). 
A minimum burn-in length of 500 steps is applied to ensure 
that burn-in is truly complete. Once this initial phase is over 
we use a further 100 steps to recalculate the parameter jump 
lengths before beginning the real Markov Chain of 10000 
accepted steps; with the acceptance rate of 25 percent rec- 



ommended for the Metropolis-Hastings algorithm ( Tegmark 
et al. 2004) this gives an effective chain length of 40000 
steps. Our set of final parameters is taken to be the me- 
dian of the Markov chain, with the lcr error bars calculated 
from the values that encompass the central 68.3 percent of 
the accepted steps. We account for limb darkening using a 



non-linear treatment based on the tables of Claret (20001, 



cussed in Anderson et al. (2011 1 



interpolating the coefficients at each step in the chain. 

The inclusion of the photometric data is an important 
point. Although we fit the RM effect to the radial velocity 
data, the transit width and depth, as well as the impact pa- 
rameter, can be determined from the photometric transit. 
These parameters have a role to play in the characterisation 
of the form of the RM anomaly. The transit width helps 
to determine the duration of the anomaly whilst the depth 
gives the planetary and stellar radii. The radii and impact 
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parameter in turn help to determine vsml, upon which the 
amplitude of the anomaly depends ( Queloz et al.|2000 1 . Al- 
though characterisation of the RM effect can be carried out 
using the spectroscopic data alone, by taking the photo- 
metric data into account in this way we ensure consistency 
across the full set of system parameters. To account for stel- 
lar jitter we initially assign a value of 1ms" 1 , below the 
level of precision of the spectrographs used for this work, 
which we added in quadrature to the in-transit photometric 
data. 

We separate our RV data by instrument, and within 
those distinctions also treat spectroscopic data taken on 
nights featuring planetary transits as separate datasets. Our 
model for the orbital RV signature treats the sets of data as 
independent, producing individual offsets and radial veloc- 
ity trends for each one. The reported solution is that for the 
set of RV data covering the greatest phase range. For com- 
pleteness, we also repeated our analysis using only RV data 
taken during nights that featured a transit event, but found 
little to distinguish them from our analysis of the the full 
set of data. 

For our RM model we use the analytic formula of |Hi-| 
rano et al. (20111. This method requires prior knowledge 



of several broadening coefficients, specifically the macro- 
turbulence, for which our estimates are noted in Tablell] 
and the Lorentzian (7) and Gaussian (/}) spectral line dis- 
persions. The line dispersions were dictated by our use of 
the HARPS instrument, which has a spectral resolution of 
R — 115000, implying an instrumental Gaussian dispersion 
of 2.61kms _1 . This was combined with the intrinsic Doppler 
linewidth, including appropriate thermal and turbulent mo- 
tion for each star, to obtain values of /3 — 3.1 kms -1 for 
WASP-16 and WASP-25, and f3 = 3.3 km s " 1 for WASP- 
31. We assumed 7 = 0.9 km s -1 in line with Hirano et al. 
and also assumed that the coefficient of differential rotation, 
a = 0. WASP-16 and WASP-25 are both slow rotators, and 
whilst WASP-31 should be considered a moderately fast ro- 
tator, without knowledge of the inclination of the stellar 
rotation axis it is difficult to place a value of a. 

We apply several Bayesian priors to \ 2 to account for 
previously known information: a prior on the eccentricity, 
allowing for the forcing of circular solutions; a prior on the 
spectroscopic vsml, using updated values of v sin I derived 
from the newly acquired HARPS spectra and the macro- 



turbulence calibration of Bruntt et al. (20101, and a prior 



enforcing a main sequence (MS) mass-radius relationship. 
This MS prior is based on that discussed in Collier Cameron 
et al. (20071, but is only applied to the stellar radius. The 



stellar mass is estimated using the calibration of |Enoch et| 

aDpoTm. 



To distinguish between models that use different combi- 
nations of priors we minimize the reduced x 2 for the spectro- 
scopic data; in cases where there is little to choose between 
the different sets of input conditions we gravitate towards 
the model with the fewest free parameters. In what follows 
we refer to % 2 as the combined x 2 for the complete data 
set, Xrv as the value for the spectroscopic RV data only, 
and Xred as the reduced x 2 f° r the spectroscopic data alone. 
Note also that we refer to the projected spin-orbit misalign- 
ment angle as A, as is more common in the literature, not ft 



4 ROSSITER-MCLAUGHLIN RESULTS 
4.1 WASP-16 



WASP-16b (Lister et al. 20091 (hereafter L09l is a close 



Jupiter analog orbiting a Solar-type star with a period of 
3.12 days. The planet is somewhat less massive than Jupiter 
but of comparable radius, whilst the host star is similar in 
mass, radius and metallicity to the Sun, but exhibits signif- 
icant lithium depletion. Our updated spectroscopic analysis 
using the HARPS spectra yields a projected stellar rotation 
velocity of vsml = 2.3 ± 0.4 km s -1 . 

Our original estimate of stellar jitter produced fits with 
1.6, leading us to re-estimate the jitter following 



X 



Wright (20051. We calculated line strengths for the calcium 



H and K emission lines in each of the HARPS spectra, and 
used these to estimate values for the chromospheric activ- 
ity metric S. These were then calibrated against the Mount 
Wilson sample (see e.g. Baliunas et al. ( 1995 1), and absolute 
magnitudes of the stars were calculated using Gray ( 1992 1 



We eventually adopted the 20th percentile value of 3.6 m s _1 
as a conservative estimate of the jitter. 

Removing the requirement for the system to obey a 



main sequence mass-radius relationship (equation 6 in Col- 
|lier Cameron et al.| ( |2007| )) produced changes of between 
and 2 percent in the stellar mass and radius, leading to in- 
creases in the stellar density of between 1 and 4 percent, for 
no discernible improvement in fit. Comparing impact pa- 
rameter values, we find that we obtain an average value of 

5+0.03 



prior active. The parameter S ( |Collier Cameron et al.|2 007), 
R, - Ro 



-2 In P(M,,i?,) 



(1) 



used to measure the discrepancy between the stellar radius 
from the (J-H) colour and that returned by the MCMC al- 
gorithm, increases from an average of 0.17 to 0.34 when the 
prior is removed, a relatively small increase as suggested 
by the modest changes in stellar parameters. We therefore 
find little to distinguish between the cases with the MS re- 
quirement applied, and those with the stellar radius freely 
varying, and choose not to apply this prior in our final so- 
lution. 

Adding a long-term, linear RV trend produced no 
improvement in Xredi an d with a magnitude of |-y < 
3ms _1 yr^ 1 we disregard the possibility that there is a such 
a trend in the spectroscopic data. Adding a prior on the 
spectroscopic vsml similarly gave almost no difference in 
the quality of the fit obtained. For most combinations of 
priors our analysis returned vsml ~ 1.2 ± 0.3kms~ , sig- 
nificantly slower than the spectroscopic value. 

Allowing the eccentricity to float again led to no signif- 
icant improvement in the fit, and all of the values of e re- 
turned by our various combinations of priors were consistent 
with e — to within 2a. We tested these small eccentricity 



values using equation 27 of Lucy & Sweeney (19711, which 



as used by T10 (strictly A = -0). 



adopts a null hypothesis of a circular orbit and considers 
an orbit to be eccentric if this is rejected at the 5 percent 
significance level. This F-test indicated that none of the ec- 
centricities were significant, and thus that a circular orbit is 
favoured. 

We therefore adopt as our conclusive solution the case 
without the MS prior active, with no prior on vsml, no 
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Orbital phase 



Figure 1. Results from the fit to the data for WASP-f6 using e = 0, no long-term radial velocity trend, no prior on the spectroscopic 
vs\nl, and without forcing the mass-radius relationship. Black, filled triangles represent data from CORALIE. Blue, filled squares 
represent data from the first HARPS run. Red, filled circles represent data from the second HARPS run. The best-fitting model is 
plotted as a solid black line. Top left: Complete radial velocity reflex motion curve. Bottom left: Residuals from theRV fit, exhibiting no 
correlation with phase. Top right: Close up of the transit region from the radial velocity curve showing the RM effect, along with the 
residuals. Bottom right: Residuals for the RV data within the RM window. 



long-term trend in velocity and a circular orbit, but we 
stress that changing the priors had little impact on the 
parameter values returned by the MCMC algorithm. Our 
adopted solution returns values of A = — 4.2° 1 +/ g 1,0 and 
vsinl = 1.2+05 km s _ > * ms i s significantly slower than the 
spectroscopic value of v sin / that we obtained from spec- 
tral analysis. However an alternative analysis of the HARPS 



spectra using the calibration of Gray (2008) provides an es- 
timate of vsinl = 1.2 ± 0.5, in good agreement with the 
value that we found from our model. Our solution also indi- 
cates a high impact parameter of 0.82+.q'q2 that reduces the 
likelihood of a degeneracy developing between A and v sin I. 
Examination of Fig. |2b| highlights this, with a triangular dis- 
tribution that is centred close to A = 0°. The main section 
of this distribution lies within the limits |A| < 20°, providing 
further evidence for the well-aligned system that was sug- 
gested by our best-fitting RM angle. From |L09] we note that 
the host star has Teff = 5700 ± 150 K, which places it in the 
'cool' category of |W10| an aligned orbit therefore fits their 
hypothesis quite nicely. 

As previously noted, the amplitude of the RM anomaly 
for WASP-16 is quite small. The aligned nature of the sys- 
tem suggests that this can be put down to the star being 
an old, slowly rotating star, which would be consistent with 



the age estimate reported by |L09| which suggests an age 
> 5 Gyr based on a lack of detectable lithium. A second 
possible explanation could be that we are in fact viewing 
the host star almost pole-on, which could still be consis- 
tent with an orbit that is aligned in the plane of the sky. 
This would lead to a low projected rotation velocity, and a 
transit across the pole of the star would have a small RM 
amplitude, as observed here. The minimum stellar inclina- 
tion is limited by the observed lithium depletion^] but such 
a structure would imply a younger age for the star owing 
to the rapi d tru e stellar rotation. Interestingly isochronal 
implies an age of 2.3+2 2 Gyr, lower than 



L09 



analysis in 

the limit implied by the lithium depletion. However new 

isochronal fits, using our results and a range of stellar mod- 



els, retu rned ages of 4.7^3 Gyr (Padova mod els; Marigo et 




), 4.8^3.3 Gyr (Yonsei-Yale mo dels; 
), 6.0+1.55 Gyr (Teramo models; 



(200"4|J^nd 5.0tH Gyr (VRSS models 



s; Pietr 
Vandei 



Dcma rque et 
infcrni et al. 



andcnBcrg, Berg- 



The abundance of lithium gives us a minimum age, as stated. If 



we assume that gyrochronology is applicable, then this provides a 



maximum true stellar rotation velocity. This in turn allows us to 
use the detected v sin / to calculate the minimum possible stellar 
inclination. 
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Table 2. A comparison of the \ 2 an d « ~ d lvalues for WASP-16 for each combination of Bayesian priors. All values of x 2 include the 
Baycsian penalties applicable for that combination of priors. 



v sin / prior 


MS 


off 


off 


off 


on 


off 


off 


off 


on 


off 


off 


off 


on 


off 


off 


off 


on 


2.3 ±0.4 


off 


2.3 ±0.4 


on 


2.3 ±0.4 


off 


2.3 ±0.4 


on 


2.3 ±0.4 


off 


2.3 ±0.4 


on 


2.3 ±0.4 


off 


2.3 ±0.4 


on 



7/ms 1 yr 



eccentricity i>sin//kms 



A/° 



X.RV 



Y 2 






0.009 


+0.010 
-0.006 


1.2 ±0.4 


9 ,+10.5 

z --n.o 





0.004 


+0.006 
-0.002 


1 2+ - 4 
l- Z -0.5 


9 n + 10.7 
~ Z - 8 -ll.l 


1 +0 ' 8 
1,u -0.8 


0.011 


+0.009 
-0.007 


1 o+ - 5 
i- u -o.6 


9 c + 13.6 
z -°-16.6 


6+ - 5 
u ' D -0.3 


0.007 


+0.007 
-0.005 


1 I -1 " 04 

i-i-o.e 


, fi +10.9 
cS - -14.8 










i 1+0.5 

i-i-o.e 


"•'-19.2 










1 2+ - 4 

1-^-0.5 


_, 9+11.0 
^■ z -13.9 


0.1 ±0.1 







1 j+0.5 


c O+10.6 
_b - 8 -14.5 
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busch & Dowlcr (2006 1). These ages further support the case 



for a slowly rotating host star, and are consistent with the 
star's observed lithium abundance. 

Careful analysis of the HARPS spectra allowed us to 
measure the chromospheric Ca II H & K emission. We find 
that \og(R' HK ) — — 5.10±0.15, indicating a low level of chro- 
mospheric activity. This rules out the possibility that the 
star is misaligned along the line-of-sight, as we would ex- 
pect much greater calcium emission from a young, rapidly 
rotating star. We note that this agrees with the work of 
|Schlaufman| (|2010|), who finds no evidence for line-of-sight 



is best suited to analysing hot, rapidly rotating exoplanet 
host stars. WASP-16 A exhibits neither of these attributes, 
but analysis is ongoing (Miller et al., in prep) and indica- 
tions are that it gives similar results for the obliquity angle 
of this system. An independent detection of the RM effect, 
also suggesting alignment, was announced at IAU Sympo- 
sium 276 by Winn, and we look forward to the published 
results with interest. 



4.2 WASP-25 



mis alignm ent in the WASP-16 system. Following |Watson et| WASP-25b ([Enoch et aLl[20TT| (hereafter [eTT| 



al. 



(20101 we calculate P ro t — 30.2_ 3 ' 8 days, which implies 



-1.2 , 



an age of 3.8_ 0g Gyr for WASP-16 according to the gy- 
rochronology method of Barnes ( 2007 1 using the updated co- 



efficients of Meibom, Mathieu & Stassun ( 2009 1 and James 
|et al.| ( |2010|>. A recent reanalysis of the WASP-1 and WASP 



2 systems ( Albrecht et al.|2011 1 highlighted the fact that in 
systems with low amplitude, low S/N RM anomalies, the an- 
gles reported tend towards 0° and 180° owing to the greater 
probability density in the distribution for A. The same study 
cautions readers against drawing strong conclusions of align- 
ment in such cases. Our data for WASP-16 certainly show 
some of the characteristics discussed in the Albrecht et al. 
study, and we have indeed found a well-aligned system with 
A close to 0. 

However there are other methods by which the align- 
ment angle of a planetary orbit can be deduced. Doppler 
tomography is an established method for mapping veloc- 
ity variations in binary stars (e.g. | Albrecht et al.| (2007, 
2009)), but its application to transiting exoplanets is in 



is a sig- 
nificantly bloated, sub- Jupiter mass planet orbiting a solar- 
type, somewhat metal-poor host star with an orbital period 
of 3.76 days. A full set of results from our analysis is dis- 
played in Table|3] One RV measurement was found to lie at 
3<7 from the best-fitting model, and to be consistent with 
the out of transit RV curve. This datum was omitted from 
our analysis, and will be discussed further later. 

We found that allowing the eccentricity to float led to a 
negligible difference in Xred> an< i that the eccentricity values 
being found were within 2<r of 0. We therefore concluded 
that the small eccentricity values being returned were aris- 
ing owing to the biases inherent in the MCMC method ( Ford 
|2006[ ), and that the orbit of WASP-25 is circular. In this we 
agree with |E11[ We confirmed this conclusion regarding a 



circular orbit using the F-test of Lucy & Sweeney (19711, 



its infancy. The technique has, to date, been used to study 
HD189733 ([Collier Cameron et al.|20J"ol), WASP-33 (|Collier| 
Cameron et al.|2010[ ) and WASP-3 ( |Miller et al.|2010[ ), arid 



which returned very high probabilities of the small eccen- 
tricity values having arisen by chance. 

We found little difference between the quality of fit for 
the equivalent cases with the MS mass-radius relation forced 
on the system, and those without the same constraint. The 
relaxation of this prior leads to larger values of A, but also in- 
creases the discrepancy between the stellar mass and radius 
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Figure 2. Posterior probability distributions derived from the Markov chains, for the fit to WASP-16 described in Fig.JT] The white 
contours mark the 62.87 percent confidence regions, the black, dashed contours the 95.45 percent confidence regions, and the black, dotted 
contours the 99.73 percent confidence regions. Marginalised, ID distributions are displayed in the side panels. Left: b and A. Right: vs'ml 
and A. This distribution has a triangular shape, and A = falls within the central body of the distribution. Both distributions have 
poorly constrained 99.73 percent confidence regions, and show a slight bias towards negative values of A. 



values. The stellar mass value varies little between runs, but 
relaxing the MS prior reduces the stellar radius by between 2 
and 3 percent, dependent on the other priors being applied. 
This leads to an increase in the stellar density of between 7 
and 12 percent from p*,MS ~ l-22p to p ti „„MS ~ 1.34/90, 
averaged across all combinations of the other priors. Con- 
sidering the impact parameter, we find that relaxing the MS 
requirement gives a value of b — 0.381q'22, whilst using the 
prior returns b = 0.44l 012 , both averaged across all other 
combinations of priors. The S parameter increases from an 
average of 3.56 to 5.92 when the prior is removed. In light 
of these differences, we elect to apply the MS prior in our 
final analysis. 

Adding a long-term linear trend in RV improved the 
Xspec of the solution, but the value of the trend varied signif- 
icantly between runs, ranging from « 2 to w 105 ms -1 yr _1 . 
We also found that in some cases the models produced when 
a trend was applied showed a notable offset from the RV 
data in transit. To check whether a trend was truly present 
in the system, 2 additional RV measurements were obtained 
using HARPS on 2010 August 25 and 26. Analysing these 
in conjunction with previously obtained data shows no ev- 
idence for a long-term RV trend, and so we disregard this 
possibility for our final solution. Introducing a prior on the 
spectroscopic v sin I produced no improvement to the qual- 
ity of fit to the data, irrespective of the other flags. We do 
not therefore apply such a prior in our final solution, and 
take this opportunity to obtain a separate measurement of 
the projected stellar rotation speed. 

Taking the results of these investigations into account, 
we select the solution with e = 0, no long-term linear trend 
in RV and no prior on v sin 7, with the main sequence mass- 



radius relation enforced. This gives A = 14.6° ±6.7, a detec- 
tion of the RM effect at 2.2<r from 0. We also obtain a value 
for the stellar rotation of vsinl = 2.9 ± 0.3 km s -1 , slightly 
greater than but in agreement with our updated spectro- 
scopic value of 2.6 ± 0.4kms _1 . The impact parameter for 
this solution is 0.44 ± 0.04. No correlation is apparent be- 
tween usinl and A (see figure[4b|, although there is evidence 
for a correlation between the impact parameter and A (see 
figure 4a |. It is possible that this correlation is responsible 



for the poor fit of the model to some parts of the RM data. 

The mechanism responsible for the outlier that we omit- 
ted from our analysis is unknown, although we note that 
|Simpson et al.|f2010[ ) experienced a similar situation in their 
analysis of the WASP-38 system, positing seeing changes 
of telescope guiding faults as possible causes. We suggest 
a third mechanism; the discrepant point might be caused 
by the planet traversing a stellar spot. In such a situation 
the spot would mask the presence of the planet, causing the 
RV measurement to diverge from the standard RM anomaly 
pattern. This scenario was suggested to explain a similar 



anomaly in the data for the WASP-3 system (Tripathi et 
|al.||20lti| , but we note that the divergence from the RM 
effect in that case showed a gradual rise and fall rather 
than the delta function change observed here, and was even- 
tually attributed to the effect of moonlight. Unfortunately 
we lack simultaneous photometry from the night of the ob- 
served spectroscopic transit, which would show the presence 
of such a spot. It is also possible that some form of transient 
event, such as a white light stellar flare, is responsible for 
the drastic, sudden change in measured RV for this point, 
although. Such events were discussed in the context of LQ 
Hya (Montes et al. 19991, and were observed to produce 
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Table 3. A comparison of the \ 2 an d X^ed values for WASP-25 for each combination of Bayesian priors. All values of x 2 include the 
Baycsian penalties applicable for that combination of priors. 



v sin I prior MS prior 7/ms x yr eccentricity i)sin//kms 1 A/ c 



~X-RV 



Y 2 

A red 



off 


off 





0011 +o.oio 

u - ui± -0.008 


off 


on 
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u.u±o_ oog 


off 


off 


103.8l 2 |j 


013+ 0014 
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off 


on 


-io.3±is : ; 


0011 +0.013 
u.un_ 008 


off 


off 








off 


on 








off 


off 


0,6 1+ 28 - 7 

aD - x -26.6 





off 


on 


2 4+ - 4 
z -^-0.3 





2.6 ±0.4 


off 





013+ 0014 
u.u±o_ oog 


2.6 ±0.4 


on 





011+ 0011 
u - uil -0 .008 


2.6 ±0.4 


off 


WOAtlH 


001 o+0.014 
u.uio_ oog 


2.6 ±0.4 


on 


Q7 1+28.0 
y '- i -25.8 
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off 
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on 
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off 


^■s + _f 5 i 





2.6 ±0.4 


on 


95.4 ±26.5 






2.8 ±0.3 
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15.9t r r i 


2.9 ±0.3 


16.8ltj 
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2.8 ±0.2 


17.li?;J 


2.8 ±0.2 


14.5±?;2 



14200 ± 169 
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14200 ± 169 
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14202 ± 169 

14185 ± 168 
14189 ± 168 
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103 ± 14 
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104 ± 14 
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91 ± 13 
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103 ± 14 

104 ± 14 

90 ± 13 

91 ± 14 
104 ± 14 
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91 ± 13 
91 ± 13 



1.3 ±0.2 
1.3 ±0.2 
1.1 ±0.2 
1.3 ±0.2 
1.3 ±0.2 
1.3 ±0.2 
1.1 ±0.2 
1.3 ±0.2 
1.3 ±0.2 
1.3 ±0.2 
1.1 ±0.2 
1.1 ±0.2 
1.3 ±0.2 
1.3 ±0.2 
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Orbital phase 







Orbital phase 



Figure 3. Results from the fit to the data for WASP-25 of our optimal solution: a circular orbit, no long-term RV trend and no prior 
on the spectroscopic vs'ml. The main sequence mass-radius relation was not enforced. The point denoted by the open square was found 
to lie 3 a from the best-fitting model, and was not included in the analysis. Legend as for Fig.JT] 
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Figure 4. Posterior probability distributions, derived from the Markov chain, for the fit to the data for WASP-25 described in Fig. [3] 
Key as for Fig.p] Left: b and A. There are appears to be some small level of correlation between the two parameters. Right: v sin I and 
A. A = falls outwith the 68.27 percent confidence contour, but within the 95.45 percent confidence contour, indicating a moderately 
significant detection of A. 



chrornospheric disturbance in the core of otherwise normal 
spectral lines. This dilution of the spectral lines could affect 
the continuum level during the flare event, and potentially 
lead to anomalous redshifting for a short period of time. 
Such an event would have to be very short duration how- 
ever, and coincide with the planet's path across the stellar 
disc. 



Should we consider WASP-25 to be aligned? |WT0] put 
forward a criterion of A Js 10° to > 3<r for misalignment; our 



result for WASP-25 clearly fails this test. |T10| suggest an 
alternative criterion of A > 30° as the limit above which we 
can be sure a system is misaligned given the average mag- 
nitude of the errors in A that are found by analysis of the 
RM effect. WASP-25 also misses this target by some mar- 
gin. But the data for the RM effect appear to be slightly 
asymmetric in Fig. |3b| suggesting that the system is mis- 
aligned (although we note that the best fitting model does 
not reflect this). 

This slight asymmetry in the RM anomaly might arise 
as a result of some form of systematic effect. We have already 
mentioned the possibility of stellar spots in the context of 
the anomalous datum omitted from our analysis. Could they 
also provide a possible explanation for the asymmetry? Con- 
sider a star on which stellar spots are more numerous in 
one hemisphere than the other during the planetary transit, 
but on which they lie away from the transit chord. As the 
planet transits the more spotty hemisphere it will hide a 
comparatively larger fraction of the photosphere, and there- 
fore mask a greater contribution to the overall flux than 
when it is transiting the less spotty hemisphere. The half 
of the anomaly corresponding to the spotted hemisphere 
would therefore have a greater amplitude than the half of 
the anomaly corresponding to the unspotted hemisphere, 



leading to an asymmetric RM effect. If the difference in the 
number and/or size of spots between the two hemispheres is 
small then the asymmetry would be only minor. This inter- 



esting systematic was discussed by Albrecht et al. (2011 1 for 



the case of WASP-2, and also seems to have played a role 



in the analysis of the CoRoT-2 RM in Bouchy et al. ( 2008 1 



In the case of WASP-25 the approaching, blue-shifted hemi- 
sphere would be required to have a slightly greater density 
of stellar spots than the receding, red-shifted hemisphere, 
which would also lead back to the possibility of a transient 
event being responsible for the anomalous datum. 

We will return to the question of WASP-25's alignment 
in section[6] 



4.3 WASP-31 



WASP-31 ( |Anderson et al.|2011[ ) is a bloated, 0.5 Mj planet 
orbiting an F-type star of sub-solar metallicity with a period 
of 3.5 days. The host star is a moderately rapid rotator, with 
wsin/ = 8.1 ± 0.5 from spectroscopy. Full results of our 
analysis can be found in Table|4] 

We found no difference between the Xred values for any 
combination of priors. We found that imposing the main- 
sequence mass-radius relation had little effect on the fit to 
the spectroscopic data, but had a deleterious effect on the 
stellar parameters. Removing the prior produced an increase 
in stellar radius of between 3 and 6 percent and a decrease 
in the stellar mass of between 1 and 2 percent, leading to a 
decrease in stellar density of between 8 and 15 percent from 
p*,MS ~ 0.67po to p t , „o ms ~ 0.62/90, averaged across all 
other combinations of priors. Comparing the impact param- 
eter and S statistic, we find b = 0.79±qq| and S = 10.2 
with no MS prior applied, and b = 0.77±o'.o4 with S = 2.9 
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Table 4. A comparison of the \ 2 an d Xred va l ues f° r WASP-31 for each combination of Bayesian priors. All values of y 2 include the 
Baycsian penalties applicable for that combination of priors. 



v sin I prior MS prior 7/ms yr 



eccentricity iisin//kms 



A/° 



~X-RV 



Y 2 



off 


off 


off 


on 


off 


off 


off 


on 


off 


off 


off 


on 


off 


off 


off 


on 


8.1 ±0.5 


off 


8.1 ±0.5 


on 


8.1 ±0.5 


off 


8.1 ±0.5 


on 


8.1 ±0.5 


off 


8.1 ±0.5 


on 


8.1 ±0.5 


off 


8.1 ±0.5 


on 
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8.1 
8.8 
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6 1+ 8 ' 1 

12 6+ 8 ' 4 





6.1 
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-0 1+ 92 
u - i -6.6 

3 4+5-7 





6 1+ 10 - 3 

5 4+ 7 ' 7 
°- 4 -8.5 



027+ - 032 
u ' uz '-0.019 

orm+ - 029 

u.uji_ 019 

023+ 0031 
u.uzo_ 017 

o 7 +0.035 
u ' uo ' -0.016 



n23+ - 029 
u.uzj_ 017 

nf)41 +0.033 

u.u4i_ 027 
0.022ig;g» 
0.038+Toil 



7.5 ±0.8 

7 7 +0.9 
'•'-0.8 

7.4 ±0.7 

7.8 ±0.8 

7.5 ±0.7 
7.5 ±0.7 
7.5 ±0.7 
7.3^ 

7.9 ±0.4 
8.0 ±0.5 

-g + 0.4 

'• 9 -0.5 

8.0 ±0.4 
7.9 ±0.4 
7.8 ±0.4 

7.8 ±0.4 

7.9 ±0.4 



2 
3.(5 



+1.1 

-2.9 

+2.9 

3.5 



j+0.7 
-0.6 



2.8+11 

o i+3.0 
°- -2.8 

2.7 ±3.0 

2.8 ±3.1 



2.1 
3.0 

2.5 
3.2 

2.7 



+2.9 

2.7 
+3.4 

3.1 
+2.8 

2.6 
+3.0 

2.9 
+2.9 

-2.7 



3.0 ±2.7 

2.; 
3.0^ 
2.7J 
3.0" 1 



d+2.7 
5 -2.9 
-,+3.0 

-2.9 
7+2.7 

-2.9 
1+3.0 

-2.9 



14703 ± 171 


64 ± 11 


0.9 ±0.2 


14708 ± 172 


64 ± 11 


0.9 ±0.2 


14700 ± 171 


63 ± 11 


0.9 ±0.2 


14695 ± 171 


63 ± 11 


0.9 ±0.2 


14702 ± 171 


64 ± 11 


0.9 ±0.2 


14706 ± 172 


64 ± 11 


0.9 ±0.2 


14698 ± 171 


63 ± 11 


0.9 ±0.2 


14698 ± 171 


64 ± 11 


0.9 ±0.2 


14693 ± 171 


64 ± 11 


0.9 ±0.2 


14703 ± 171 


64 ± 11 


0.9 ±0.2 


14698 ± 171 


63 ± 11 


0.9 ±0.2 


14702 ± 171 


64 ± 11 


0.9 ±0.2 


14697 ± 171 


64 ± 11 


0.9 ±0.2 


14701 ± 171 


65 ± 11 


0.9 ±0.2 


14701 ± 171 


64 ± 11 


0.9 ±0.2 


14705 ± 171 


65 ± 11 


0.9 ±0.2 



when the requirement for the star to be on the MS is en- 
forced. Owing to the much more favourable 5* statistic, and 
the influence on the stellar parameters, we elect to use re- 
sults which account for the MS relationship. Adding a linear 
velocity trend gave no discernible difference in the quality 
of the fit to the spectroscopic data, and with a magnitude 
of 1 7 1 < 13ms _1 yr _1 we conclude that no such trend is 
present in the system. Adding a prior on the spectroscopic 
v sin I made little difference to the results despite the rela- 
tively rapid rotation, so we again choose the simpler route 
and neglect such a prior. Finally, we choose a circular solu- 



tion; the F-test of Lucy & Sweeney ( 1971 1 shows that the 



small eccentricity values returned when e is allowed to float 
are insignificant. 

Our optimal solution is therefore that obtained with 
no vsinl prior, no velocity trend, the MS prior active, and 
e = 0. This set of priors gives A = 2.8° ± 3.1 , leading to 
the conclusion that the WASP-31 system is well-aligned. It 
is worth noting that this would be the conclusion whichever 
combination of priors we chose, as all of the values of A that 
we obtained lie within 1.2a of 0°. The impact parameter is 
0.771q'o2- The stellar rotation for this solution has a value 
of v sin I — 7.5±0.7kms _1 . As with our result for WASP-16 
this is slower than the spectroscopic value, but in this case 
the value agrees to within lcr. Again, an alternative analy- 



sis using the calibration of Gray ( 1992 I returns a value of 
vsinl (7.5 ± 0.5 km s -1 ) more similar to our MCMC result. 



WASP-31 is not included in the sample of Schlaufman ( 2010 1 



owing to its time of publication. In order to check the possi- 
bility of misalignment along the line-of-sight, we follow the 
method of Schlaufman and calculate the rotation statistic, 
O. The age of WASP-31 A is somewhat uncertain however; 
its lithium abundance, gyrochronology and the presence of 
a close companion all suggest ages of ~ 1 Gyr, whilst pre- 



vious stellar model fits imply an older age of 4 ± 1 Gyr. 
We reassess the isochronal fit for the system, obtaining 
ages of 4.0ligGyr (Padova models), 2.8tJ; Gyr (Yonsei- 
Yale models), 3.5l 13 Gyr (Teramo models) and 2. 8+12 Gyr 
(VRSS models). Using these estimates we calculate values 
of 9 = -4.5, -3.0, -3.7 and -2.8 respectively; WASP-31 is 
therefore rotating more rapidly than expected given its age 
in both cases. The chance of significant misalignment along 
the line-of-sight therefore seems slim; the inclination of the 
WASP-31 b's orbit is 84.6 ± 0.2°, leaving little room for an 
increase in rotation velocity owing to line-of-sight misalign- 
ment. 



INTEGRATION INTO THE ENSEMBLE OF 
RESULTS 



The analysis of |W10' provides a good starting point for in- 
tegrating our new results into the existing ensemble of RM 
measurements. Fig.[7]reproduces fig. 2 from their paper, with 
the addition of all complete RM measurements made since 
its publication (except WASP-23 (jTriaud et "aL"1|2011[), for 



which the result is still highly uncertain, and WASP-26 JAn 



derson et al.|2011), which showed only a very low amplitude 



and was classed as a non-detection); we list these planets in 
Table[5] We also elect to include most of the systems that 
|W10| disregard during their analysis as having insufficiently 
precise measurements of Arlin order to provide a full picture 



2 HAT-P-2, C0R0T-I, CoRoT-3, HD149026, Kepler-8, TrES-1 
and TrES-2. See references within [WlOl Although WASP-2 has a 
measured value for A, the most recent analysis of the system failed 
to detect a signal I jAlbrecht et al.|2011[ | and thus we continue to 
exclude this system. 
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Figure 5. Results from the fit to the data for our adopted solution for WASP-31, with a circular orbit, no prior on the spectroscopic 
vsml, no long-term radial velocity trend, and the mass-radius relationship applied. Legend as for Fig.ll] 
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Figure 6. Posterior probability distributions, derived from the Markov chain, for the fit to the data for WASP-31 described in Fig. pj] 
Key as for Fig. 'reffig:W16prob. Left: b and A. Right: vsinl and A. A = lies well within the main body of the distribution. 
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Figure 7. Projected stellar obliquity, A, as a function of stel- 
lar effective temperature for all systems with confirmed measure- 
ments. WASP-16 is marked by a green, filled triangle, WASP-25 
by a blue, filled circle, and WASP-31 by a red, filled square. The 
vertical dotted line marks the distinction between 'cool' and 'hot' 
systems, whilst the horizontal dotted line marks 30°, the angle 
above which a system is considered to be misaligned in |W10| and 

PTTol 



of the current state of RM analysis. Whilst it is true that 
making a definitive statement regarding alignment is more 
difficult for these systems owing to their large uncertainties, 
the criteria for granting misaligned status should take ac- 
count of this. We are also interested in comparing our new 
measurements to the general form of the current ensemble. 
Omitting the systems listed above does not simplify this 
task, so we elect to include them. 

WASP-31 has an effective temperature of 6300 ± 100 K, 
which falls with lcr of the border between the 'hot' and 'cool' 
categories of |W10| albeit tending towards the 'hot' side. We 
cannot therefore draw any conclusions as to how it affects 
the trend proposed in that paper. 

With an effective temperature of 5750 ± 100 K, WASP- 
25 falls into the 'cool' category (Teff ^ 6250 K) of |W10 
which they find to be preferentially aligned 



their sample 

gives a probability of misalignment for 'cool' stars of 0.17. 
Updating this result using our expanded sample changes the 
probability to 0.20 using either the criterion ofjWIO, or to 
0.13 using the criterion of |A| > 30° from |T10| It is worth 
noting here that the apparently large differences in misalign- 
ment probability between the two criteria are an artefact of 
the sample size, which is still relatively small at 48 systems 
(30 'cool', 18 'hot'). Switching between the two criteria only 
changes the number of aligned systems by two for the 'cool' 
sub-sample, and has no effect on the number of misaligned 
systems in the 'hot' sub-sample. Under both criteria the ap- 
parent alignment of WASP-25 b's orbit is in accordance with 
the [WTO] hypothesis. WASP-16, T eff = 5700 ± 150 is also 
classified as a 'cool' system. All available information points 
towards this system being well-aligned, and it therefore fits 
well with the hypothesis of |W10| 

The final interesting point about Fig. [7] is the apparent 
lack of systems with mildly retrograde, close to polar orbits. 
There are currently only two systems with 80° ^ A ^ 110°, 
and only one more with 110° ^ A ^ 140° This relatively un- 



populated region is less noticeable when considering ip owing 
to the increased size of the error bars, but it is still apparent. 
We speculate that truly polar orbits are perhaps unstable for 
some reason. Or perhaps it is simply our inability to deter- 
mine the inclination of the stellar rotation axis that is at 
fault. It may be that some 'aligned' systems actually have 
close to polar orbits if this angle is accounted for. It may also 
be that we simply have yet to observe very many systems in 
this region of the parameter space, and future publications 
may provide the data to fill this underpopulated area. 

It has not been remarked upon before in this context, 
but a drop in the number of systems at mid-range obliq- 
uity angles is clearly predicted by the theoretical ip angular 



distribution of Fabrycky & Tremaine (2007). It also clearly 



shows up in the angular distribution for the complete set 
of known obliquity angles, fig. 10 in |T10| We reproduce this 
figure in Fig.lSJ adding the probability distributions of the 
planets in TablelS] as well as those of the planetary systems 
from this study, ip, the true misalignment angle, is given by 



cos tp = cos I cos i + sin I sin i cos A, 



(2) 



where I is the inclination of the stellar rotation axis to the 
line-of-sight, and i is the inclination of the orbital axis to 
the line-of-sight. To calculate the V distribution for each 
planet we carried out 10 6 Monte Carlo simulations, drawing 
values for / from a uniform cos/ distribution to represent 
the case in which stellar rotation axes are randomly ori- 
ented on the sky. We also accounted for the error bars on i 
and A by drawing values from a Gaussian distribution with 
our optimal solution values as the mean values, and scaled 
to the uncertainties in those values. The individual plan- 
ets' distributions were then summed to produce our total 
distribution, which is similar to that of |T10| and still com- 
pares favourably to the theoretical histogram from |Fabrycky| 
& Tremaine (20071. The drop in probability at mid-range 



angles is in line with the underpopulated region of Fig.: 7: 
and our additions bring the primary, low-angle peak closer 
in shape to the theoretical distribution. The overall shape 
of the secondary peak is less clear; it is still dominated by 
contributions from individual systems owing to the smaller 
number of planets with strongly misaligned orbits as com- 
pared to the number of aligned or weakly misaligned sys- 
tems, but appears as though it may be broader and more 
shallow than the theoretical prediction. 

Fig. [8] requires the assumption that the /, the stellar 
inclination, is isotropic and that the angular distribution is 
unimodal. However the discussion of |W10| implies that the 
distribution is in fact bimodal. A clearer demonstration of 
the agreement between theoretical predictions and current 
observations is therefore to look at the distribution in A. 
This requires the converse transformation of the predicted 



t/j distribution of Fabrycky & Tremaine ( 2007 1 into A 



We reproduce the lower panel of fig. 9 from |T10"l taking 
into account the additional measurements of A from Tabled 
For HAT-P-7 and HAT-P-14, both of which have published 
A > 180.0°, we used the negative angle equivalent (360— A). 
This cumulative A distribution avoids both of the assump- 
tions inherent in Fig.[8] Agreement between the observa- 
tional data and the theoretical predictions of Fabrycky <fe 
(20071 has been improved, particularly for low- to 



Tremaine 



mid-range angles, but the observational data are still slightly 
lacking in high obliquity systems compared to the theoret- 



The alignment of WASP- 16, WASP-25 and WASP-31 13 

Table 5. Relevant data for the planetary systems for which the Rossiter-McLaughlin effect has been characterised since the publication 
of |W10| We add these systems to the |Winn et al.| sample to bring the ensemble of results up to date and allow us to better analyse the 
place of WASP-25 and WASP-31 within that ensemble. 



System 



*/° 



usiniykms" 1 Teff/K 



A/° 



Reference 



CoRoT-18 

HAT-P-4 

HAT-P-6 

HAT-P-8 

HAT-P-9 

HAT-P-11 

HAT-P-14 

HAT-P-16 

HAT-P-23 

HAT-P-30 

KOI-13.01 

WASP-1 

WASP-7 

WASP-19 

WASP-22 

WASP-24 

WASP-38 

XO-3 

XO-4 



86.5+JI 
88.76±S;H 

85.51 ±0.35 
87.5 ± 1.4 

86.5 ±0.2 
89.17±g:« 

83.52 ±0.22 

86.6 ±0.7 
85.1 ± 1.5 
83.6 ±0.4 
85.0 ±0.4 
90 ±2 
87.2±?:§ 
79, 4 ±0.4 
88.26 ±0.91 
83.64 ±0.29 
88.83^ 
82.5 ± 1.5 
88.8 ±0.6 



j+0.51 
-0.55 



-,+0.95 
-0.56 



8.0 ± 1.0 

5.83 ±0.35 
7.5 ± 1.6 
14.5 ±0.8 
12.5 ± 1.8 
1.00J 
8.18 ±0.49 
3.9 ±0.8 

7.8 ± 1.6 
3.07 ±0.24 

65 ± 10 

7 +i.4 

U -'-0.5 
14 ±2 

4.63 ±0.26 
4.42 ± 0.34 
7.0 ±0.9 
8.58 ±0.39 
18.4 ±0.2 

8.9 ±0.5 



5440 ± 100 
5860 ± 80 
6570 ± 80 
6200 ± 80 
6350 ± 150 
4780 ± 50 
6600 ± 90 
6158 ± 80 
5905 ± 80 
6304 ± 88 
8511 ±400 
6110 ±45 
6400 ± 100 
5500 ± 100 
5958 ± 98 
6075 ± 100 
6150 ± 80 
6430 ± 50 
6397 ± 70 



3+26 
-10 



10 ±20 
4.9 ±11.9 
166 ± 10 

_ 17 +9.2 

-16 ±8 

103^ 

189.1 ±5.1 

-10.0 ±16 

15 ±22 

73.5 ±9.0 

23 ±4 

-59^ 

86 ±6 

4.6 ±5.2 

22 ±16 

-4.7 ±4.0 

15+ 33 
10 43 

37.4 ±2.2 
-46.7 ±7.1 



1+99 
-26 
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Figure 8. The total distribution of the true obliquity angle, if>, 
for the complete sample of systems for which A has been mea- 
sured. The dotted histogram represents the theoretical distribu- 
tion of |Fabrycky fc Tremaine| ( |2007fr . The dashed line represents 
the limit of ij> = 30° above which a system is considered to be 
misaligned. The overall forms are comparable, and the total ip 
distribution is similar to fig. 10 of |T10| The shape of the primary 
peak agrees well with theoretical predictions. The overall shape of 
the secondary, high angle peak in the distribution is less clear, but 
may be more shallow and broader than anticipated. The sudden 
drop in probability density at mid-range angles, around ip fa 90°, 
has become more pronounced when compared to the distribution 

of rnoi 




Figure 9. Cumulative probability histogram for A. The solid line 
denotes observational data, whilst the dotted line denotes the the- 
oretical distribution of |Fabrycky &; Tremaine| p007), converted 
from ip to A. The vertical blue, dashed line marks A = 30°, the 
limit above which a planetary orbit is considered to be misaligned. 
The agreement between the two distributions has improved with 
the addition of measurements made since the publication of |T10| 
particularly at mid-range angle, prograde orbits, but the obser- 
vational data is still lacking in high obliquity systems compared 
to the theoretical prediction. 



ical histogram, whilst showing more low-obliquity systems 
than expected. 
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6 A NEW MISALIGNMENT TEST 

The Rossiter-McLaughlin effect has now been measured for 
48 transiting exoplanets, but as of yet there seems little con- 
sensus as to the best way of classifying them as aligned or 
misaligned. For most of the systems with measurements of 
A this is not a serious problem; either |A| > 90°, or the error 
bars are such that the obliquity is consistent with zero. But 
as the number of RM measurements continues to grow, there 
will be an increasing number of systems in a similar situa- 
tion to WASP-25, which exhibits a mildly asymmetrical RM 
anomaly but does not fulfil any of the current misalignment 
criteria. 

There are two main criteria currently in use by the com- 



munity. W10 use |A| > 0° at ^ 3a significance to define a 



misaligned system. T10 take |A| > 30° as their threshold, 



on the basis that errors in the obliquity angle are of the or- 
der of 10°, and therefore this gives 3<r significance as well. 
We would like to introduce a new test for misalignment that 
takes a completely different approach to these. 

We consider the set of WASP planets for which the RM 
effect has been characterised using RV data, including the 
systems presented in this study. We neglect the WASP-33 
system for which the misalignment angle has been measured 



only through Doppler tomography (Collier Cameron et al 



2010), and disregard the ambiguous results for WASP-23 



(Triaud et al. 20111 and WASP-2 (Albrecht et al. 2011) 



For reasons of consistency we use the RV based solution of 



Tripathi et al. (2010) for our intial conditions for WASP-3, 



rather than the more recent tomographical study of Miller et 
al. ( 2010 1. The full set of planets sample is listed in Table|6] 
Our test is based on the Bayesian Information Criterion 
(BIC) ( |Liddle|2007| ), 



BIC = x!v + fcln(n), 



(3) 



where k is the number of parameters and n is the number 
of data. Changing the value of A only affects the form of 
the model RV curve in-transit; we therefore just consider 
those RV points that lie within a region of the RV curve 
around phase defined by the fractional transit width when 
computing the second term of the BIC. The number of pa- 
rameters changes according to the choice of priors applied 
to the MCMC run; adding a long-term RV trend, fitting the 
RM effect, and allowing the eccentricity to float all add one 
or more additional parameters to the model. 

We carry out two MCMC analyses for each of the sys- 
tems in our sample, using the same combination of priors 
for both. The first analysis allows both \/v sin I cos A and 
^/v sin / sin A to float, whilst the second forces an aligned 
orbit by fixing ^/v sin I sin A = 0. We calculate the BIC for 
both runs, before calculating B — BIC a jj„ n /BIC. For the 3 
systems presented herein we use our adopted solutions, and 
carry out an additional run to provide the aligned case. We 
plot the results for all of the systems as a function of the 
sky-projected alignment angle. 

We find several distinct groups of systems within our 
results, which lead us to define three categories of align- 
ment into which systems with RM measurements can be 
classified. Five systems, including WASP-16 and WASP-31, 
were found to haveB ^ 0.980, implying that the model with 
A = provides a better fit that the free-floating A model. Of 
these five systems, all would be classified as aligned accord- 



ing to either of the existing misalignment criteria. A fur- 
ther four systems, including WASP-25, are clustered around 
B — 1.00, forming a distinct group in figure fTOb] Forcing 
an aligned orbit would seem to make little difference to the 
quality of the fit between data and model in these cases. 
Of these systems three would clearly be classed as aligned 
according to|W10|and|T10| but the fourth (WASP-1) would 



actually be classed as misaligned according to |W10| The 
remaining systems clearly lie distinct from those discussed 
so far, and many are clearly classifiable as misaligned, with 
|A| > 100° and B > 1.5. 

In light of these results, we define three categories of 
alignment. Systems for which B ^ 0.99 we classify as mis- 
aligned. Those with B ^ 1.01 we classify as aligned. Systems 
falling between these categories, with 0.99 < B < 1.01 we 
classify as of indeterminate alignment. We would also de- 
fine a fourth category, that of 'no detection', as containing 
those systems with vsinl consistent with to within la, 
but our current sample contains no systems that meet this 
requirement. 

Some of the systems in TableRj] warrant a little more 
examination. WASP-16, despite the relatively poor quality 
of the RM data that we obtained, can be more strongly con- 
sidered aligned than WASP-31 with its high quality data. 
This is an interesting, if puzzling result, but does provide 
further evidence to support our previous conclusion of an 
aligned system for WASP-16. WASP-25 is classified as un- 
determined under our new criteria, possibly owing to the 
relatively poor match between the shape of the RM anomaly 
and the best-fitting model. However we note that it lies very 
close to the boundary between the 'undetermined' and 'mis- 
aligned' categories. Our new MCMC runs for WASP-1 and 
WASP-4 produce very large error bars on A, but end up in 
different categories despite both failing the |T10| criterion of 
misalignment. Examining their respective RM anomalies we 
note that both have very low amplitudes, but that the data 
for WASP-4 is of significantly better quality than that for 
WASP-1. It is likely that this is responsible for the differ- 
ence in classification. In addition, |T10| noted a substantial 
correlation between A and usin/ for WASP-4, arising due 
to the low impact parameter, which may be producing the 
large lower error. WASP-38 also exhibits a significant error 
bar on A, and we again note that the RM data is again of 
somewhat poor quality. New observations of WASP-38 using 
HARPS may help to improve the quality of the results for 
the system, allowing us to draw firmer conclusions (Brown 
et al, in prep.). 



7 CONCLUSIONS 

We have presented analysis of the Rossiter-McLaughlin ef- 
fect for WASP-16, WASP-25 and WASP-31. We find WASP- 
16 to have a very low amplitude signal, but the use of two 
complete spectroscopic transits has enabled us to determine 
a sky-projected alignment angle of A = — 4.2°J~ 1 3 1 g . For 
WASP-25 we find a mildly asymmetric RM anomaly with 
A = 14.6° ± 6.7, and for WASP-31 we obtain A = 2.8° ± 3.1, 
indicating a well-aligned system. 

Since WASP-31 lies so close to the effective temperature 
that divides the classes of 'hot' and 'cool' planet hosts, we 
cannot say how its alignment affects the pattern proposed 
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Table 6. Relevant data for our new misalignment criterion, for a sample of WASP planets with existing Rossiter-McLaughlin measure- 
ments. A values are those obtained from our new MCMC analyses. BIC values were calculated from the spectroscopic x 2 values, using the 
number of in-transit RV measurements only. Our new misalignment criterion defines systems with a BIC ratio B ^ 1.01 as misaligned, 
those with B ^ 0.99 as aligned, and those with 0.99 < B < 1.01 as of indeterminate status. 



System 



reference 



A/° 



v sin //km s 



BIC 



BIC 



align 



ABIC B 



WASP-1 
WASP-3 


Albrecht et al.|(J2011 
frIplitSTet = ai = [(|2TM 


) 60.2±i 6 3 6 
37.9±»i 3 8 

42 0+ 14 - 3 

26.2+J;* 

7 c+20.9 
'•°-19.1 

) 85-Ot 9 8 ;o 
-106.71^ 

-28.0+J;° 

-lss.si^ 7 

-134.5+ ? - 3 
20-5±i?;|' 

-i-«±S:S 

-6.9^: 4 8 

-6.ll 3 8 3 7 


1.3 ±0.5 
12 n+l.l 

9 c+0.4 

3.5 ±0.2 

i 7+0.3 

J -'-o.a 

26.3lij 

9 S+ 04 
A8 -0.3 

2.8 ±0.3 

4 5+0-4 
4 -°-0.3 

9.8 ±0.3 
12.9 ±0.3 
3.2 ±0.2 

5 1+0-4 
°- i -0.3 

8.2 ±0.3 


255.2 ±22.6 
294.7 ±24.3 
86.8 ± 13.2 

186.0 ±19.3 

134.7 ±16.4 

285.8 ±23.9 

380.3 ±27.8 

171.2 ±18.5 

154.4 ±17.6 

342.7 ±26.2 

118.8 ±15.4 
81.5 ± 12.8 

123.1 ±15.7 

241.3 ±22.0 


256.7 ±22.7 
308.0 ±24.8 
91.3 ± 13.5 
199.2 ±20.0 
132.0 ±16.2 
451.7 ±30.1 
1092.5 ±46.7 
193.5 ±19.7 

555.7 ±33.3 
986.9 ±44.4 
119.0 ±15.4 
79.8 ± 12.6 

119.8 ±15.5 
240.8 ±21.9 


1.5 

13.3 

4.5 

12.1 

-2.7 

165.9 

712.2 

22.3 

401.3 

644.2 

0.2 

-1.7 

-3.3 

-0.5 


1.006 

1.045 


WASP-4 


Triaud et al. 


2010 


) 
) 


1.052 


WASP-5 


Triaud et al. 


2010 


1.071 


WASP-6 


Gillon et al. ( 


2009 


0.980 


WASP-7 


Albrecht et a 


^20 


12 


1.580 


WASP-8 


Queloz et al. 


2010 


) 

) 
) 
) 


2.873 


WASP- 14 


Joshi et al. 1 1 


0091 


1.130 


WASP-15 


Triaud et al. 


2010 


3.599 


WASP-17 


Triaud et al. 


2010 


2.880 


WASP-18 


Triaud et al. 


2010 


1.002 


WASP-19 


Hellier et al. 


2011 


0.979 


WASP-24 


Simpson et a' 


. 120: 





0.973 


WASP-38 


Simpson et a 


. (2010 


0.998 



WASP-16 this study 
WASP-25 this study 
WASP-31 this study 



-4.2+ 
14.6 ± 
2.8 ± J 



Ll.O 

13.9 


1 2+ - 4 


115.6 ±15.2 


112.1 ±15.0 


-3.5 


0.970 


6.7 


2.9 ±0.3 


116.5 ±15.3 


117.5 ±15.3 


0.8 


1.009 


.1 


7.5 ±0.7 


73.7 ± 12.1 


72.2 ± 12.0 


-1.5 


0.980 



by |W10| WASP-25 on the other hand at first appears to 
strengthen their hypothesis, with the existing misalignment 
criteria of both WIO] and |T10| labelling it aligned. We have 
also presented a new method for determining the alignment 
or otherwise of an exoplanetary orbit. Our test is based on 
the BIC statistic, and bases the misalignment or alignment 
of a system on the ratio of the values of the BIC for the free A 
case and the aligned case. We classify systems with B ^ 1.01 
as misaligned, those with B ^ 0.99 as aligned, and those 
with 0.99 ^ B ^ 1.01 as of indeterminate classification. 
WASP-25 falls in this last category, albeit very close to the 
boundary with the 'misaligned' classification. 

The results presented herein bring the analysis of the 
ensemble of systems with confirmed stellar obliquities up to 
date. Our results have done little to change the overall pic- 
ture presented by |T10| instead strengthening the agreement 
with theoretical predictions for the distributions of both the 
projected and true stellar obliquities. We should not be too 
hasty to assume that we have solved the problem of hot 
Jupiter migration however; new discoveries are constantly 
causing us to re-evaluate our current understanding. 
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Figure 10. B as a function of A for the sample of planets in TableRj] as well as the systems presented in this study. WASP-16 is denoted 
by a filled, green triangle. WASP-25 is denoted by a filled, blue circle. WASP-31 is denoted by a filled, red square. The horizontal dotted 
line marks B = 1.00. The two vertical dotted lines denote |A| = 0° and |A| = 30°, the existing criterion for misalignment. Left: All 
data. Right: A close-up of the heavily populated region in the lower left of the plot, around B = 1.00 and |A| = 0°. This shows the 
separation of the systems into several distinct groupings, which lead us to define three categories of alignment. This changes the existing 
classification of some systems. 
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Table Bl. Radial velocity data for WASP-16 obtained using the 
CORALIE high precision cchcllc spectrograph. 



HJD(-2450000) 


Tcxp/s 


RV/kms- 1 


cr RV /km s 1 


4535.864842 


900 


-1.99772 


0.01591 


4537.849158 


1427 


-1.96688 


0.00853 


4538.858364 


1800 


-2.00734 


0.00899 


4558.780835 


1800 


-1.83336 


0.00723 


4560.709473 


1800 


-2.00513 


0.00725 


4561.688137 


1800 


-1.82730 


0.00785 


4589.705102 


1800 


-1.84255 


0.00875 


4591.706755 


1800 


-2.03571 


0.00892 


4652.495906 


1800 


-1.82493 


0.00808 


4656.551645 


1800 


-2.02421 


0.00787 


4657.577293 


1800 


-1.96640 


0.00957 


4663.539741 


1800 


-2.02961 


0.00969 


4664.616769 


1800 


-1.78590 


0.01108 


4682.521501 


1800 


-1.98118 


0.00754 


4881.869213 


1800 


-2.02245 


0.00813 


4882.801025 


1800 


-1.83289 


0.00823 


4884.737094 


1800 


-2.04565 


0.00778 


4891.805707 


1800 


-1.90043 


0.00798 


4892.723980 


1800 


-1.83413 


0.00891 


4941.728231 


1800 


-1.88737 


0.00748 


4943.730102 


1800 


-2.04677 


0.00753 


4944.739293 


1800 


-1.91359 


0.00860 


4945.799895 


1800 


-1.85815 


0.00807 


4947.745317 


1800 


-1.93960 


0.00741 


4948.673112 


1800 


-1.82992 


0.00743 


4972.707323 


1800 


-1.93123 


0.00854 


4975.733486 


1800 


-1.93144 


0.01100 


4982.647535 


1800 


-1.83433 


0.01036 


4984.642389 


1800 


-2.04210 


0.00892 


4985.694776 


1800 


-1.81561 


0.00802 


5391.544362 


1800 


-1.80313 


0.00867 
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Table Al. Parameters from the best-fitting, adopted models for the three WASP planetary systems studied in the main text 



System v sin/ prior MS prior 7/ms 1 yr 1 eccentricity jitter /ms 1 t)sin//kms 1 X/° 



WASP-16 


off 


off 








1.0 


1 2+ - 4 


_ 4 0+11.0 


WASP-25 


off 


on 








3.6 


2.9 ±0.3 


14.6 ±6.7 


WASP-31 


off 


on 








1.0 


7.5 ±0.7 


2.8 ±3.1 



Table B2. Radial velocity data for WASP-16, for the first transit 
obtained using the HARPS high precision echcllc spectrograph on 
the night of 2010 March 21. 



Table B3. Radial velocity data for WASP-16, for the second 
transit obtained using the HARPS high precision cchelle spectro- 
graph on the night of 2011 May 12. 



HJD(-2450000) 


Tcxp/s 


RV/kms" 1 


°"ijv/ kms 1 


HJD(-2450000) 


Tcxp/s 


RV/kms- 1 


°\Rv/kms J 


5275.661171 


1800 


-1.80610 


0.00337 


5685.845943 


900 


-2.02724 


0.00305 


5275.907691 


1800 


-1.78144 


0.00188 


5687.838150 


900 


-1.79259 


0.00383 


5276.661941 


500 


-1.88264 


0.00533 


5692.662149 


900 


-1.99841 


0.00380 


5276.668446 


500 


-1.89263 


0.00546 


5692.796210 


900 


-1.96847 


0.00309 


5276.674824 


500 


-1.88914 


0.00583 


5693.517817 


900 


-1.81013 


0.00298 


5276.681375 


500 


-1.87845 


0.00547 


5693.800775 


900 


-1.78196 


0.00285 


5276.687753 


500 


-1.89900 


0.00557 


5694.581176 


600 


-1.88349 


0.00344 


5276.694258 


500 


-1.88947 


0.00569 


5694.588340 


600 


-1.88597 


0.00302 


5276.700693 


500 


-1.89868 


0.00555 


5694.595389 


600 


-1.88670 


0.00310 


5276.707094 


500 


-1.89750 


0.00552 


5694.602900 


600 


-1.88871 


0.00305 


5276.713599 


500 


-1.88945 


0.00586 


5694.610180 


600 


-1.89619 


0.00290 


5276.720046 


500 


-1.91121 


0.00614 


5694.616904 


500 


-1.89308 


0.00323 


5276.726493 


500 


-1.89650 


0.00623 


5694.623386 


500 


-1.88547 


0.00323 


5276.732929 


500 


-1.88385 


0.00657 


5694.629531 


500 


-1.89159 


0.00309 


5276.739376 


500 


-1.90596 


0.00640 


5694.635631 


500 


-1.89935 


0.00312 


5276.745812 


500 


-1.90686 


0.00704 


5694.641904 


500 


-1.89371 


0.00307 


5276.752143 


500 


-1.90101 


0.00698 


5694.648003 


500 


-1.89694 


0.00318 


5276.758579 


500 


-1.91736 


0.00742 


5694.654149 


500 


-1.91191 


0.00298 


5276.765605 


500 


-1.90950 


0.00627 


5694.660364 


500 


-1.91005 


0.00311 


5276.771589 


500 


-1.91143 


0.00447 


5694.666406 


500 


-1.91009 


0.00300 


5276.778140 


500 


-1.91692 


0.00439 


5694.672609 


500 


-1.91343 


0.00311 


5276.784344 


500 


-1.91573 


0.00479 


5694.678824 


500 


-1.91547 


0.00327 


5276.790838 


500 


-1.92779 


0.00526 


5694.684924 


500 


-1.91384 


0.00313 


5276.797459 


500 


-1.92243 


0.00510 


5694.691070 


500 


-1.91484 


0.00330 


5276.803964 


500 


-1.90902 


0.00433 


5694.697227 


500 


-1.91701 


0.00308 


5276.810411 


500 


-1.92567 


0.00391 


5694.703373 


500 


-1.91926 


0.00308 


5276.816441 


500 


-1.91540 


0.00416 


5694.709460 


500 


-1.91211 


0.00335 


5276.823178 


500 


-1.92366 


0.00424 


5694.715664 


500 


-1.91534 


0.00353 


5276.829336 


500 


-1.92742 


0.00454 


5694.721821 


500 


-1.91948 


0.00349 


5276.835887 


500 


-1.92693 


0.00522 


5694.727979 


500 


-1.92566 


0.00346 


5276.842334 


500 


-1.92962 


0.00572 


5694.734078 


500 


-1.92109 


0.00367 


5276.848723 


500 


-1.94183 


0.00699 


5694.740351 


500 


-1.92435 


0.00377 


5276.855228 


500 


-1.94543 


0.00926 


5694.746612 


500 


-1.91926 


0.00364 


5276.861907 


500 


-1.92508 


0.00830 


5694.752596 


500 


-1.93024 


0.00354 


5277.630948 


1800 


-2.02847 

"i r\r\ or yi 


0.00222 


5695.501446 


900 


-2.03120 


0.00281 


5277.861599 


1800 


— 1.99854 


0.00196 










5278.632376 


1800 


-1.82733 


0.00398 










5278.857922 


1800 


-1.79546 


0.00208 










5279.627285 


1800 


-1.84379 


0.00264 










5279.913540 


1500 


-1.91554 


0.00242 










5280.624797 


1800 


-2.03079 


0.00266 










5280.916481 


1200 


-2.00824 


0.00283 
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Table B4. Radial velocity data for WASP-25 obtained using the 
CORALIE high precision echclle spectrograph. 



Table B5. Radial velocity data for WASP-25 obtained using the 
HARPS high precision echelle spectrograph. The point denoted 
by * was omitted from the analysis (see text for details). 



HJD(-2450000) tcxp/s RV/kms" 



°ftv/kms 



HJD(-2450000) tcxp/s RV/kms" 1 a RV /kms -1 



4829.822664 


1800 


-2.57717 


0.01282 










4896.769798 


1800 


-2.65105 


0.01069 


5296.540546 


1200 


-2.54661 


0.00329 


4940.709168 


1800 


-2.71589 


0.01154 


5296.635060 


1200 


-2.54381 


0.00398 


4941.704336 


1800 


-2.61855 


0.01153 


5297.506446 


1200 


-2.61464 


0.00509 


4942.725717 


1800 


-2.57632 


0.01238 


5297.518749 


400 


-2.62531 


0.01007 


4943.637434 


1800 


-2.61828 


0.01246 


5297.523714 


400 


-2.63325 


0.00942 


4944.715466 


1800 


-2.67966 


0.01207 


5297.528714 


400 


-2.60898 


0.01012 


4945.726530 


1800 


-2.61467 


0.01305 


5297.533807 


400 


-2.61031 


0.01041 


4946.616622 


1800 


-2.58169 


0.01266 


5297.538714 


400 


-2.60809 


0.00999 


4947.601618 


1800 


-2.64133 


0.01096 


5297.543714 


100 


-2.59973 


0.01056 


4947.791245 


1800 


-2.68927 


0.01347 


5297.548668 


400 


-2.59327 


0.01108 


4948.613002 


1800 


-2.70418 


0.01098 


5297.553761 


400 


-2.58984 


0.01028 


4949.803142 


560 


-2.55132 


0.01819 


5297.559131 


400 


-2.61140 


0.01786 


4950.622083 


1800 


-2.59141 


0.01348 


5297.563761 


400 


-2.60108 


0.01177 


4951.695324 


1800 


-2.70149 


0.01218 


5297.568668 


400 


-2.60927 


0.01087 


4971.645302 


1800 


-2.67821 


0.02101 


5297.573715 


400 


-2.60539 


0.01075 


4972.672436 


1800 


-2.56129 


0.01319 


5297.578761 


400 


-2.62964 


0.01153 


4973.515713 


1800 


-2.58676 


0.01269 


5297.583668 


400 


-2.62703 


0.01118 


4974.678659 


1800 


-2.71413 


0.01359 


5297.588669 


400 


-2.61584 


0.01225 


4975.537940 


1800 


-2.66695 


0.01384 


5297.593715 


400 


-2.64141 


0.01190 


4976.683662 


1800 


-2.55567 


0.01304 


5297.598669 


400 


-2.64658 


0.01232 


4982.619435 


1800 


-2.66448 


0.02096 


5297.603715 


400 


-2.65755 


0.01234 


4983.621314 


1800 


-2.56777 


0.01486 


5297.608715 


400 


-2.67520 


0.01246 


4983.644577 


1800 


-2.59698 


0.01450 


5297.613761 


400 


-2.67558 


0.01254 


4984.578450 


1800 


-2.55837 


0.01474 


5297.618715 


400 


-2.68567 


0.01244 


4985.609967 


1800 


-2.69905 


0.01189 


5297.623669* 


400 


-2.63635 


0.01215 


4995.555496 


1800 


-2.50858 


0.01396 


5297.628854 


400 


-2.67450 


0.01065 


5009.628712 


1800 


-2.60564 


0.01823 


5297.633761 


400 


-2.65389 


0.00837 


5010.596729 


1800 


-2.53871 


0.02313 


5297.638715 


400 


-2.63022 


0.00840 










5297.643773 


400 


—2.63126 


0.00885 










5297.648727 


400 


-2.61768 


0.00871 










5297.653727 


400 


-2.63157 


0.00862 










5297.658681 


400 


-2.63982 


0.00841 










5297.663773 


400 


-2.62371 


0.00834 










5297.668773 


400 


-2.63776 


0.00800 










5297.673727 


400 


-2.64716 


0.00811 










5297.678727 


400 


-2.63753 


0.00780 










5297.683681 


400 


-2.63915 


0.00781 










5297.688820 


400 


-2.63818 


0.00823 










5297.693727 


400 


-2.64295 


0.00763 










5297.698774 


400 


-2.62870 


0.00772 










5297.703635 


100 


-2.63667 


0.00757 










5297.708727 


400 


-2.63667 


0.00781 










5297.713727 


400 


-2.63031 


0.00803 










5297.718681 


400 


-2.65168 


0.00756 










5297.723774 


400 


-2.64533 


0.00785 










5297.833578 


1200 


-2.65676 


0.00352 










5298.535157 


1200 


-2.69608 


0.00406 










5298.716015 


1200 


-2.69119 


0.00285 










5298.830796 


1200 


-2.69107 


0.00287 










5299.544943 


1200 


-2.60603 


0.00327 










5299.701761 


1200 


-2.62922 


0.01842 










5299.838220 


1384 


-2.57573 


0.01224 
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Table B6. Radial velocity data for WASP-31 obtained using the 
CORALIE high precision echclle spectrograph. 



Table B7. Radial velocity data for WASP-31 obtained using the 
HARPS high precision echelle spectrograph. 



HJD(-2450000) 


texp/s 


RV/kms- 1 


o-flv/kms ! 


HJD(-2450000) 


texp/s 


RV/kms- 1 


o-jjv/kms 1 


4835.809755 


1800 


-0.20260 


0.02945 


5298.496133 


1200 


-0.15638 


0.00896 


4837.773728 


1800 


-0.08457 


0.03163 


5298.749145 


1200 


-0.17218 


0.00931 


4840.765776 


1800 


-0.07974 


0.03752 


5299.504441 


1200 


-0.17141 


0.01114 


4880.767231 


1800 


-0.20651 


0.03722 


5299.716991 


1200 


-0.29541 


0.04800 


4939.627676 


1800 


-0.07528 


0.03471 


5300.509357 


1200 


-0.07910 


0.01138 


4941.567460 


1800 


-0.19372 


0.02721 


5300.742948 


1200 


-0.08587 


0.01161 


4942.654757 


1800 


-0.12324 


0.02598 


5301.582822 


1200 


-0.11516 


0.00934 


4943.610624 


1800 


-0.07939 


0.03435 


5301.597544 


1200 


-0.11250 


0.01003 


4944.555415 


1800 


-0.16578 


0.02620 


5301.612960 


1200 


-0.12527 


0.00984 


4945.544750 


1800 


-0.16158 


0.03220 


5301.627126 


900 


-0.10841 


0.01418 


4946.591778 


1800 


-0.04422 


0.03441 


5301.638028 


900 


-0.09394 


0.01141 


4947.555094 


1800 


-0.14075 


0.02656 


5301.648803 


900 


-0.06309 


0.01227 


4948.588069 


1800 


-0.18274 


0.02861 


5301.659393 


900 


-0.08776 


0.01100 


4950.597039 


1800 


-0.08667 


0.03009 


5301.670307 


900 


-0.11321 


0.01136 


4951.608264 


1800 


-0.22821 


0.03366 


5301.681082 


900 


-0.13542 


0.01119 


4971.548671 


1800 


-0.02806 


0.07763 


5301.691973 


900 


-0.16717 


0.01218 


4973.489951 


1800 


-0.09189 


0.03084 


5301.702354 


900 


-0.19579 


0.01297 


4974.608541 


1800 


-0.02130 


0.03225 


5301.713349 


900 


-0.16796 


0.01407 


4975.511060 


1800 


-0.12628 


0.03381 


5301.723522 


900 


-0.16142 


0.01884 


4983.595539 


1800 


-0.10028 


0.04415 


5301.734934 


900 


-0.14695 


0.01771 


4984.467744 


1800 


-0.06112 


0.02895 


5301.750477 


1200 


-0.13378 


0.01244 


4985.530406 


1800 


-0.18886 


0.03032 


5301.764794 


1200 


-0.14088 


0.01586 


4994.508045 


1800 


-0.09735 


0.03237 


5301.780360 


1200 


-0.15031 


0.01738 


4994.531307 


1800 


-0.09484 


0.03469 


5305.613529 


1200 


-0.18166 


0.01096 


4995.463376 
4995.486741 


1800 


-0.12246 


0.03547 


5307.584640 


1200 


-0.06475 


0.00809 


1800 


—0.20999 


0.03214 










4996.459605 


1800 


-0.19034 


0.03472 










4996.482971 


1800 


-0.10984 


0.03117 










4999.536757 


1800 


-0.21980 


0.05522 










4999.560099 


1800 


-0.15670 


0.06560 










5006.521354 


1800 


-0.18493 


0.04217 










5012.492297 


1800 


-0.09667 


0.03739 










5013.497045 


1800 


-0.19082 


0.04459 










5029.465780 


1800 


-0.11660 


0.05041 










5168.846768 


1800 


-0.15389 


0.01852 










5203.782854 


2700 


-0.20152 


0.02084 










5290.715577 


2700 


-0.09351 


0.01809 










5291.699859 


2700 


-0.12273 


0.01821 










5293.696819 


2700 


-0.07455 


0.01987 










5294.733907 


2700 


-0.09785 


0.01925 










5296.704921 


2700 


-0.13168 


0.01772 










5298.693406 


2700 


-0.18210 


0.01718 










5300.589643 


2700 


-0.06706 


0.02250 










5326.628560 


2700 


-0.13831 


0.01945 










5327.604475 


2700 


-0.08093 


0.02737 










5328.608442 


2700 


-0.09912 


0.01913 










5334.544675 


2700 


-0.05571 


0.02129 











